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Abstract. 

For a lattice A with n vertices and dimension d equal or higher than two, the number 

of spanning trees A^st(A) grows asymptotically as exp(nzA) in the thermodynamic 

C^ , limit. We present exact integral expressions for the asymptotic growth constant za 

S ' for spanning trees on several lattices. By taking different unit cells in the calculation, 

' ■ many integration identities can be obtained. We also give 2:a(p) on the homeomorphic 

(-H I expansion of fc-regular lattices with p vertices inserted on each edge. 
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1. Introduction 

A problem of fundamental interest in mathematics and physics is the enumeration of 
the number of spanning trees on the graph G, Nst{G). This number can be calculated 
in several ways, including as a determinant of the Laplacian matrix of G and as a special 
case of the Tutte polynomial of G P . A previous study on the enumeration of spanning 
trees and the calculation of their asymptotic growth constants was carried out in Ref. 
|2j. In that work, closed- form integrals for these quantities were given. In this paper 
we shall present exact integrals for the asymptotic growth constant for spanning trees 
on five Archimedean lattices and the diamond structure. The hypercubic lattices with 
both nearest-neighbor and next-nearest-neighbor edges will also be considered. We shall 
show how integration identities can be obtained with different choices of unit cells in the 
calculation. We shall also derive the asymptotic growth constant for the homeomorphic 
expansion of a regular lattice. 

2. Background and Method 

We briefly recall some definitions and background on spanning trees and the calculation 
method that we use jllini- Let G = iy,E) denote a connected graph (without loops) 
with vertex (site) and edge (bond) sets V and E. Let n = v{G) = \V\ be the number 
of vertices and e{G) = \E\ the number of edges in G. A spanning subgraph G' is a 
subgraph of G with v{G') = \V\, and a tree is a connected subgraph with no circuits. A 
spanning tree is a spanning subgraph of G that is a tree and hence e{G') = n — 1. The 
degree or coordination number ki of a vertex Vi & V is the number of edges attached to 
it. Two vertices are adjacent if they are connected by an edge. The adjacency matrix 
A{G) of G is the n x n matrix with elements A{G)ij = 1 if f j and Vj are adjacent 
and zero otherwise. The Laplacian matrix Q{G) is the n x n matrix with element 
Q{G)ij = ki5ij — A{G)ij. One of the eigenvalues of Q{G) is always zero; let us denote 
the rest as \{G)i, 1 < i < n — 1. A basic theorem is that Nst{G) = (l/n) YliZi ^{G)i 
pp. For (i-dimensional lattices A with rf > 2 in the thermodynamic limit, Nst{-^) 
grows exponentially with n as n —* oo; that is, there exists a constant z^ such that 
A^5t(A) ~ exp(n2;A) as n — i> oo. The constant describing this exponential growth is thus 
given by 

za = lim n"MnA^5r(A) (1) 

where A, when used as a subscript in this manner, implicitly refers to the thermodynamic 
limit of the lattice A. If each pairs of neighboring vertices are connected by p-iold 
edges instead of just one edge, the number of spanning trees should be multiplied by 
the factor p^~^ and the resultant asymptotic growth constant increases by the factor 
In p. Henceforth, we only consider simple graphs without multiple edges. A regular d- 
dimensional lattice is comprised of repeated unit cells, each containing u vertices. Define 
a{n,n') as the u x u matrix describing the adjacency of the vertices of the unit cells n 
and n', the elements of which are given by a(n, fi')ij = 1 if f j G n is adjacent to Vj G n 
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and zero otherwise. Assuming that a given lattice has periodic boundary conditions, 
and using the resultant translational symmetry, we have a{h,h') = a{n — n'), and we 
can therefore write a{n) = a(ni, ■ ■ ■ , hd). In Ref. |2I a method was derived to calculate 
A^5t(A) and z\ in terms of a matrix M which is determined by these a{n). For a d- 
dimensional lattice which is fc-regular that each of its vertices has the same degree k, 
define 

M(^i,---,0,) = A;-l-^a(n)e™-® (2) 

h 

where in this equation 1 is the unit matrix and stands for the (i-dimensional vector 
(^i,---,^,). Thenia 

where D{9i,- ■ ■ ,9d) = det(M(6'i, ■ ■ ■ ,9d)) is the determinant of the matrix. 

An Archimedean lattice is a uniform tiling of the plane by regular polygons in 
which all vertices are equivalent [U IH] . Such a lattice can be defined by the ordered 
sequence of polygons that one traverses in making a complete circuit around the local 
neighborhood of any vertex. This is indicated by the notation A = (HiPD' meaning 
that in this circuit, the regular pj-sided polygon occurs contiguously a^ times. All eleven 
Archimedean tilings are fc-regular. 

The number of spanning trees is the same for a planar graph G and its dual G*, 
and the number of the vertices of G* is given by the Euler relation n* = \E\ — n + 1. 
Because \E\/n = k/2 for a fc-regular graph Gk, the asymptotic growth constants of Gk 
and Gl satisfy [2j 

Consequently, zc^ = 2zg^, zqi = zq^, zqi = 2zgJ3 and zc^ = zgJ2. 

For a fc-regular graph Gk, a general upper bound is zq^ <\nk. A stronger upper 
bound for a fc-regular graph Gk with k > 3 can be obtained from the bound [HI Ej 

NsT{Gk)<(^)ibkr (5) 



where 



_^nk\nk 
(k-l)''-^ 



(6) 



[k{k-2p-^ 
With eq. (0), this then yields 

ZG, < Hh) ■ (7) 

It is of interest to see how close the exact results are to these upper bounds. For this 
purpose, we define the ratio 
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Two graph G and G' are homeomorphic to each other if one of them, say G', can 
be obtained from the other, G, by successive insertions of degree-2 vertices on edges of 
G |n|. This process is called homeomorphic expansion. We shall denote a lattice A with 
p vertices inserted on each edge as A{p). 

3. Asymptotic Growth Constants 

The asymptotic growth constants za for many d > 2 lattices have been considered by 
several authors [3111131101 HI]- For certain two-dimensional lattices, z\ can be expressed 
in terms of known functions. For instances, Zsq = 4C/7r, Ztri = (ln3)/2+(6/7r)Tz2(l/v^) 
[H 1131, and 2(4.8.8) = C/n + (ln(y2 - l))/2 + (Ti2(3 + 2\/2))/7r [IJ, where G is the 
Catalan constant and Ti2 is the inverse tangent integral function ^31. While Zhc = ztri/^, 
due to the duality between honeycomb and triangular lattices, it is non-trivial to have the 
relations Zkag = (-2trJ + ln6)/3 and ^(3.12.12) = (-2iri + ln(15))/6 given in p. In this section, 
we study the following five Archimedean lattices, (3^.4^), (3^.4.3.4), (3^.6), (3.4.6.4) and 
(4.6.12), that have not been considered before. Most of the calculations done so far 
have been performed for lattices with nearest neighbor edges. The consideration can 
be extended to lattices with edges connecting farther vertices. We shall take the d- 
dimensional hypercubic lattice with both nearest-neighbor and next-nearest-neighbor 
edges for illustration. 

3.1. (33.42) Lattice 

The (3^.4^) lattice can be constructed by starting with the square lattice and adding a 
diagonal edge connecting the vertices in, say, the upper left to the lower right corners 
of each square in every other row as shown in Fig. ^ (a). The simplest unit cell for this 
lattice contains two vertices 2/(33,42) = 2 as labeled in Fig. [H (a). We have 

M (f) f))-( 5-2cos^i _i _ e^^i - e^^^ A 

The determinant can be calculated to be 

1^(33.42) (^1,^2) = 22 - 22COS01 +4cos2^i - 2 cos ^2 - 2cos(0i - 62) (10) 
such that 

1 r de, r de^ , r^ .„ „ . 

^(33.42) = 2 / -^ / ^In P(3«.42)(^l,^2) 



2 7_. 27r 



ll-llcos0 + 2cos^^ 



+ v/(l -cos0)(7-2cos0)(17- 13 cos ^ + 2 cos2 6 

= 1.406925832 (11) 

where one integration is carried out followed by a numerical evaluation of the remaining 
integral. 
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Figure 1. (a) The (33.42) lattice, (b) The (3^4.3.4) lattice, (c) The (3''.6) lattice, (d) 
The (3.4.6.4) lattice, (e) The (4.6.12) lattice. Vertices within a unit cell are labeled. 
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3.2. (32.4.3.4) Lattice 

The (3^.4.3.4) lattice is shown in Fig. Q (b) which has the structure of a net of squares. 
Taking each square as a unit cell with z^(32 .4.3,4) = 4, we have 

^(32. 4.3.4) (^1, 6^2) 



/ 5 -1 - e*^i 

-1 - e-*^i 5 

.1 _ £-^^2 _g-«S2 

The determinant can be calculated to be 



1 - e*^2 \ 

J02 



v 



4_gi02 _g 

5 -1 - e-*^i 

-1-e^^i 5 



(12) 



^(32.4.3.4) (^1,^2) = 336 - 144(cos0i + 00862) +4:{cos'^9i + cos^ ^2) 



56 cos 6*1 cos 6^2 



(13) 



such that 



2^(32 .4.3.4) 



1 r d9i r d92 



27r 



27r 



-0(32 .4.3. 4) (6^1, 6^2) 



in2+ir^r^in 

4 ./^ 27r ./ ^ 27r 



+ 



dOi r dOi 



27r 



27r 



In 



6 + 4cos(^)cos(^) 



COS ( — ) 



COS ( — ) 



A 



92. 



6 — 4 cos( — ) cos( — ) — COS ( ^ ) — COS , 



1.410855646 . 



(14) 



3.3. (3^6) Lattice 

The (3^.6) lattice is shown in Fig. Q (c) which has the structure of a net of hexagons. 
Taking each hexagon as a unit cell with 1^(34.0) = 6, we have 

M(34.6)(ei,^2) 



/ 



v 



5 




-1 


_giei 


_^i{0i+O2) 


_gi(ei- 


-1 




5 


-1 


_gj(ei+e2) 


_gie 


_g-*ei 




-1 


5 


-1 


-e^« 


g-i(ei+e2) 


— e 


-i{ei+e2) 


-1 


5 


-1 


g-i(ei+e2) 


- 


_g-i02 


-e~*^2 


-1 


5 


-1 




-e-^^2 


_giei 


_gi6li 


-1 



316*2 



— e 



(15) 



and the determinant -0(34.6) (6*1, 6^2) is given by /(5032, 1624, 56), where we define the 
function 



/(a, b,c) = a — b[cos 61 + cos 62 + cos(6'i — ^2)] 

— c[cos 61 cos 6*2 + (cos 61 + cos 62) cos(^i 
+ 8 cos 6*1 cos 62 cos{6i — 62) . 



d2)] 



(16) 
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Thus, 



Z(3iS) 



1 r dOi r de2 



6 



27r 



^In 



27r 



/(5032,1624,56) 



1.392023563 . 



:i7) 



This result apphes to both enantiomorphic forms of the (3'^. 6) lattice. 



3.4. (3.4.6.4) Lattice 

The (3.4.6.4) lattice is shown in Fig. ^ (d) which also has the structure of a net of 
hexagons. Taking each hexagon as a unit cell with 1^(3,4.6.4) = 6, we have 

M(3.4.6.4)(^l,^2) 
/ 



The determinant D 
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-1 
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_^^e, 


-1 


4 


(3.4.6.4) (6'l, 6^2) 


is given by /( 


1144,376,8) such that 





1 r rf^i r de2 , 



/(1144,376,8) 



1.144801124 . 



(18) 



(19) 



5.5. (4.6.12) Lattice 
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The determinant D 



2^(4.6.12) 



(4.6.12) (^1,^2) is given by /(13432,4344, 136) such that 

1 r dOi r de2 



12 }_^ 27r J_^ 27r 



In 



/(13432, 4344, 136) 



0.7777955061 



(21) 



3. 6. d- Dimensional Hypercubic Lattice with Next-Nearest-Neighbor Edges 

The d-dimensional hypercubic lattice, denoted as c{d), with nearest-neighbor edges has 
been studied before in Refs. j3 HHI- Now consider the hypercubic lattice with next- 
nearest-neighbor edges included, denoted as c{d),nnn. The one-dimensional circuit 
graph with next-nearest-neighbor edges corresponds to a strip of the triangular lattice 
with width two. We have 



Mc(l),nnn{0) 



4 - 2 COS 



-1 



je 



-1 



'id 



4-2cos^ 



(22) 



The determinant -Dc(i),rmn(^) is equal to 2(1 — cos0)(7 — 2 cos 6*) such that 

Zcii),nnn = (1/2) ln[(7 + 3v^)/2] = 0.9624236501 (23) 

which agrees with Refs. [21 E]. For d > 2, the coordination number is given by 



k 



c(d),nnn 



2d + 4(2) = 2d'^. We have the general expression 



Zc{d),nnn 



f [fl— In('2rf2-2^cos0i-2^cos^icos0 



for rf > 2 . 



(24) 



The square lattice with next-nearest-neighbor edges corresponds to c{d = 2), nnn that 

"" dOi r de2 



^sq,nnn 



^ dO 



2ti 



In 



8 — 2 cos 9i — 2 cos 62 — ^ cos 9i cos 62 



In 



TT 



4 - cos^ + v/3(l-cos^)(5 + cose) 



(25) 



An exact closed form expression for this integral can be derived. We begin by recasting 
the integral in the equivalent form. 



^sq,nnn 



ln3 + - 

TX 



7r/2 



I In 



sin0 + 



(1/3) sin^ 



In 3 + 4/(1/^3) 



where 



lia) 



TT 



7r/2 



In ( sin 6 + \l 1 — a^ sin^ 



(26) 



(27) 



In eq. (jTTjl . with no loss of generality, we take a to be nonnegative. We will give a 
general result for I [a) with < a < 1 and then specialize to our case a = 1/^/S. First, 
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and 



to obtain 



we note that /(I) = C/n, where C is the Catalan constant. Next, taking the derivative 
with respect to a and doing the integral over (p in eq. ()27|) . we get 

a/2 - (2/7r) tanh"^ a 

^ ^") = iiT^) • ^^^> 

To calculate I{a), we then use I (a) — /(O) = j^ I'{x)dx and observe that 

/(O) = - / d(j) ln(sin0 + 1) = —^ + — . (29) 

TT Jq 2 71 

We also make use of the integrals 

Jt^^^ = ^M1 + «') (30) 

/■" tanh^^ X , , , _i C tt, /l + a\ 1^ /l + a\, , 

/ r— ax = tan atanh a-\ h— In — Tio (31) 

Jq 1 + x2 2 8 \l-aJ 2 \l-aJ^ ' 

CI, /(l + a^)(l-a)\ 2 T , 1 1^. /l + a\,„^s 

Ua) = — +- In ^ -^^- — ^ ) tan"^ a tanh"^ a+-Ti2 ) (32) 

^ ' vr 4 V 4(l + a) / vr vr Vl-a/^ ^ 

where Ti2(x) is the inverse tangent integral [Tlj . 

Ti2(x) = / dy 

Jo y 

= x[ ,F,{[1, 1/2, 1/2], [3/2, 3/2], -x^) ] . (33) 

Evaluating our result (|H^ for /(a) at a = 1/v^ and substituting into eq. (HHj), we 
obtain the exact, closed-form expression 

Zs<,,nnn = — + ln(2 - Vs) - - tauh-^ ^ + -Ti2(2 + V3 ) 
vr 3 V3 TT 

= 1.943739360 . (34) 

We note that Zgq and ztri for the square and triangular lattices, respectively, can be 
expressed in terms of the Ti2 function U21E1- The numerical evaluations for 3 < c? < 5 
are given by 

Zci3),nnn = 2.841787536 

^c(4),nnn = 3.442108 

Zci5),nnn = 3.898251 . (35) 

With the observation that Zc(d),nnn ~^ ln(2(i^) for d ^ oo and similar behavior for other 
(i-dimensional lattices [2], we conjecture that for a d-dimensional lattice A with degree 
kd, zj^ approaches to \Ti{kd) from below in the large d limit. 

4. Integral Identities 

As the method reviewed in section |21 is general, there is a freedom to specify the unit 
cell for the calculation. Many integral identities can be obtained with different choices 
of unit cells in the calculation of the spanning trees. 
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Consider the body- centered cubic (bcc) lattice with coordination number k^cc = 8 as the 
first example. Although the conventional bcc cell contains two vertices, the primitive cell 
contains only one vertex. In terms of the unit vectors x, y, z of the Cartesian coordinate, 
the primitive translation vectors of the bcc lattice are {x + y — z)/2, {x — y + z)/2 and 
{—X + y + z)/2. Hence, 



Zbc 



31n2 + 



27r J_^ 271 



2tx 



In 



hcc\ 



h, 02,0-3 



(36) 



where 



i^6cc(^l,^2,^3) = l-^(cOS^l+COSe2 + COS03 + COs(01+02 + ^3)) • (37) 



Comparing Zbcc with the expressions given in eq. (5.2.7) of Ref. and in eq.(16) of Ref. 
[TT] , we find 

"" dOi r dOi r dOs 

^ J_^ 27 J_^ 27 



Fbcc{dl, 02,03) 



1 °° 1 



5F4([1, 1,3/2, 3/2, 3/2], [2,2,2,2],!) 
24 



(3^ 



4-2. fee Lattice 



Next, consider the face-centered-cubic (fee) lattice with coordination number kfcc = 12. 
The conventional fee cell contains four vertices, but the primitive cell again contains only 
one vertex. The primitive translation vectors of the fee lattice are [x + y)/2, {y + z)/2 
and {z + x)/2. Hence, 



1 (,o^M r^^i r^^2 r^^3, 

.,. = ln(12)+/ ^ 1^ 1^1- 



'fcc[ 



^1,^2,^3 



(39) 



where 



FfcciOi, 62,63) = 1 - ^f cos 6^1 + cos 6^2 + cos 6*3 



+ cos( 



e. 



C0S(^2 - ^3 



cos(e3 - ^i)) . (40) 



The equality of this triple integral and the triple integral given in eqs.(18) and (19) of 
Ref. U5 is therefore established. 



4-3. Diamond Structure 

The diamond structure can be viewed as two fee structures displaced from each other by 
one-quarter of a cubic diagonal. Each vertex has four nearest-neighbor, kdiamond = 4, at 
the corners of a tetrahedron. The conventional unit cube contains eight vertices, so that 
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^diamond — 8, located at (0, 0, 0), (2, 25 0), (2, 0, 2), (0, 2, n)) (45 I) 4)5 (4) 45 4)) (45 I) 4) 



and 



i,|,|). We have 
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63) \ 



The determinant of this matrix is exactly the same as that for the fee lattice, i.e., 
12^F(^i,^2,^3) where ^(^1,^2,^3) is given in eq.(19) of Ref. [Uj. Therefore, 

Zdiamond = Zfj2 ~ 1.206459950 . (42) 

This result can be verified by taking the primitive basis with two vertices located at 



(0, 0, 0) and 



•1 1 l^ 
-4' 4' 4> 



The corresponding matrix is given by 



M, 



4 



g-j^i _ g-«6»2 _ g-«^3 



diamond 



piSi pi82 gi 



(43) 



with determinant 12Fjcc(6'i, 6'2, 6'3) where FjcdOi, 62,03) is given above in eq. (jiUI) . 



(41) 



4-4- Square Lattice 

Let us take the square lattice as another illustration. It has been pointed out in Ref. 
|2] that d = 2 case of the generalized body-centered-cubic lattice bcc{d) is the square 
lattice. The unit cell for bcc{2) with i/ = 2 is shown in Fig. El (a). If one takes the unit 
cell containing two vertices as shown in Fig|21(b), the corresponding matrix is 

4-2cosei -l-e^^2 
_l_g-ie2 4_2cos0i 



M. 



sq,2 



(44) 



with determinant 

DsqA^i, O2) = 14-16 cos ^1+4 cos^ ^i - 2 cos 62 . 



(45) 



A possible choice of the unit cell containing three vertices for the square lattice is shown 
in Fig. 12(c). The corresponding matrix is 



M, 



sg,3 



/4-2cos^i -1 

-1 4-2cos^i 



102 



-e 
-1 



\ 



V 



-ie2 



-1 



(46) 



4 - 2 cos 6^1 j 



Spanning Trees on Lattices and Integration Identities 12 

with determinant 

Dsq,3{0i,92) = 52-9Ocos0i+48cos^^i -Scos^^i -2 cos ^2 • (47) 

A possible choice of the unit cell containing four vertices for the square lattice is shown 
in Fig. 121(d). The corresponding matrix is 

/ 4 -1 - e~*^2 



M. 



sq,A 







jei 



^itix 







V 



-1 -e 




-itix 



4 
-1 - e-'^-' 



-1 -e' 


-1 - e'^'' 
4 



\ 



/ 



with determinant 



^59,4(^1,^2) = 128-64(cos^i + cos^2)+4(cos^i -cos ^2) 



In 
In 
In 



-Dsq,2(^l,^2) 



D 



sq,3l 



D 



sq,i\ 



h,02 



h,02 



(4^ 



(49) 



All of these considerations give the same Zgq for the square lattice |S1 E] , 
_AC _1 r dOi r dOi 

_ 1 r (^ r d02 

"31,271,27 

_ 1 r dOi r de2 
" 4 1, 27 y_, 27 

The same procedure can be carried on for lattices with different choices of unit cells, 
and many integral identities can be obtained. 



(50) 



5. Homeomorphic Expansion 

The method reviewed in section |21 can be applied to lattices which are not fc- regular. We 
find that for a fc-regular lattice A, the asymptotic growth constant of its homeomorphic 
expansion A(p) is given by 

;f-l)ln(p+l) + ZA 



ZK(p) 



fp+1 



(51) 



This can be understood as follows. As the number of occupied edges by a spanning tree 
on the original lattice is n — 1, the number of unoccupied edges on the original lattice 
is e{G) — {n — \) = {k/2)n — n + 1. Now on each unoccupied edge we insert p degree-2 
vertices. There are p + 1 ways to make these vertices connected in order to construct 
spanning trees on the homeomorphic expanded lattices, so that 

Nst{I^{p)) = {p + l)t"-"+iAr5r(A) . (52) 

Eq. (j51|) follows by noticing that the total number of vertices of the homeomorphic 
expanded lattice is {k/2)np + n. For the homeomorphic expansion of the square, 
triangular and honeycomb lattices, 

\n{p + 1) + Zsg 



Zsqip) 



2p + l 
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(a) 



(b) 



3 2 

4 1 



c) (d) 

Figure 2. Different unit cells for the square lattice, (a) The choice for the bcc(2) 
lattice with ly = 2. (b) The other choice with v = 2. (c) A choice with ly — 3. (d) A 
choice with v = A. Vertices within a unit cell are labeled and edges within each unit 
cell are connected by thick lines. 



^tri{p) 



^hc{p) 



2 ln(p + 1) + Ztri 

3p+l 
I ln(p + 1) + Zhc 

|p+l 



(53) 



where p is a positive integer. For a simple illustration, consider the square lattice with 
one vertex inserted on each edge. As shown in Fig. El the number of vertices in each 
unit cell i^sq(i) is equal to three. We have 



/ 



M< 



9(1) 



-1 



-26*2 



V 







-1 - e^^2 

4 
-1 - e'^' 







(54) 



The determinant Dsq{i) is equal to 8 — 4 cos 6*1 — 4 cos 6*2, which is twice of the 
corresponding expression for the regular square lattice [2,, such that 2;sg(i) = (ln2 + 

Zsq)/^- 
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Figure 3. sg(l), the homeomorphic expansion of the square lattice with p = 1 vertex 
inserted on each edge. Vertices within a unit cell are labeled and edges within each 
unit cell are connected by thick lines. 

Table 1. Values of za and ta- 



A 


d 


k 


ZA 


rA 


(33.42) 


2 


5 


1.406925832 


0.9486371781 


(32.4.3.4) 


2 


5 


1.410855646 


0.9512869040 


(34.6) 


2 


5 


1.392023563 


0.9385891391 


(3.4.6.4) 


2 


4 


1.144801124 


0.9411423250 


(4.6.12) 


2 


3 


0.7777955061 


0.9292789200 


c{l),nnn 


1 


4 


0.9624236501 


0.7912095935 


c{2),nnn 


2 


8 


1.943739360 


0.9681095375 


c{3),nnn 


3 


18 


2.841787536 


0.9933044903 


c{4:),nnn 


4 


32 


3.442108 


0.9978270 


c{5),nnn 


5 


50 


3.898251 


0.9990858 


diamond 


3 


4 


1.206459950 


0.9918321170 



6. Discussion 

We summarize the values of z^ and ta for various lattices A in Tabled Our results agree 
with the observations made in Ref. [21 . Namely, za is relatively large for large value of 
k] while for a fixed k value, z\ increases with the spatial dimension of the lattice. An 
example is given by the circuit graph with next-nearest-neighbor edges, (3.4.6.4) lattice 
and the diamond structure (all with k = A). The ratio r^ increases with dimension and 
approaches to one. 
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